Coarse-Grained Modeling of the SARS-CoV-2 Spike Glycoprotein by Physics-Informed Machine Learning

Coarse-grained (CG) modeling has defined a well-established approach to accessing greater space and time scales inaccessible to the computationally expensive all-atomic (AA) molecular dynamics (MD) simulations. Popular methods of CG follow a bottom-up architecture to match properties of fine-grained or experimental data whose development is a daunting challenge for requiring the derivation of a new set of parameters in potential calculation. We proposed a novel physics-informed machine learning (PIML) framework for a CG model and applied it, as a verification, for modeling the SARS-CoV-2 spike glycoprotein. The PIML in the proposed framework employs a force-matching scheme with which we determined the force-field parameters. Our PIML framework defines its trainable parameters as the CG force-field parameters and predicts the instantaneous forces on each CG bead, learning the force field parameters to best match the predicted forces with the reference forces. Using the learned interaction parameters, CGMD validation simulations reach the microsecond time scale with stability, at a simulation speed 40,000 times faster than the conventional AAMD. Compared with the traditional iterative approach, our framework matches the AA reference structure with better accuracy. The improved efficiency enhances the timeliness of research and development in producing long-term simulations of SARS-CoV-2 and opens avenues to help illuminate protein mechanisms and predict its environmental changes.

Coarse-Grained Modeling of the SARS-CoV-2 Spike Glycoprotein by Physics-Informed Machine Learning

Coarse-grained (CG) modeling has defined a well-established approach to accessing greater space and time scales inaccessible to the computationally expensive all-atomic (AA) molecular dynamics (MD) simulations. Popular methods of CG follow a bottom-up architecture to match properties of fine-grained or experimental data whose development is a daunting challenge for requiring the derivation of a new set of parameters in potential calculation. We proposed a novel physics-informed machine learning (PIML) framework for a CG model and applied it, as a verification, for modeling the SARS-CoV-2 spike glycoprotein. The PIML in the proposed framework employs a force-matching scheme with which we determined the force-field parameters. Our PIML framework defines its trainable parameters as the CG force-field parameters and predicts the instantaneous forces on each CG bead, learning the force field parameters to best match the predicted forces with the reference forces. Using the learned interaction parameters, CGMD validation simulations reach the microsecond time scale with stability, at a simulation speed 40,000 times faster than the conventional AAMD. Compared with the traditional iterative approach, our framework matches the AA reference structure with better accuracy. The improved efficiency enhances the timeliness of research and development in producing long-term simulations of SARS-CoV-2 and opens avenues to help illuminate protein mechanisms and predict its environmental changes.

Modeling of the thermal properties of SARS-CoV-2 S-protein.

We calculate the thermal and conformational states of the spike glycoprotein (S-protein) of SARS-CoV-2 at seven temperatures ranging from 3°C to 95°C by all-atom molecular dynamics (MD) µs-scale simulations with the objectives to understand the structural variations on the temperatures and to determine the potential phase transition while trying to correlate such findings of the S-protein with the observed properties of the SARS-CoV2. Our simulations revealed the following thermal properties of the S-protein: 1) It is structurally stable at 3°C, agreeing with observations that the virus stays active for more than two weeks in the cold supply chain; 2) Its structure varies more significantly at temperature values of 60°C-80°C; 3) The sharpest structural variations occur near 60°C, signaling a plausible critical temperature nearby; 4) The maximum deviation of the receptor-binding domain at 37°C, corroborating the anecdotal observations that the virus is most infective at 37°C; 5) The in silico data agree with reported experiments of the SARS-CoV-2 survival times from weeks to seconds by our clustering approach analysis. Our MD simulations at µs scales demonstrated the S-protein's thermodynamics of the critical states at around 60°C, and the stable and denatured states for temperatures below and above this value, respectively.

Modeling of the thermal properties of SARS-CoV-2 S-protein

We calculate the thermal and conformational states of the spike glycoprotein (S-protein) of SARS-CoV-2 at seven temperatures ranging from 3°C to 95°C by all-atom molecular dynamics (MD) µs-scale simulations with the objectives to understand the structural variations on the temperatures and to determine the potential phase transition while trying to correlate such findings of the S-protein with the observed properties of the SARS-CoV2. Our simulations revealed the following thermal properties of the S-protein: 1) It is structurally stable at 3°C, agreeing with observations that the virus stays active for more than two weeks in the cold supply chain;2) Its structure varies more significantly at temperature values of 60°C–80°C;3) The sharpest structural variations occur near 60°C, signaling a plausible critical temperature nearby;4) The maximum deviation of the receptor-binding domain at 37°C, corroborating the anecdotal observations that the virus is most infective at 37°C;5) The in silico data agree with reported experiments of the SARS-CoV-2 survival times from weeks to seconds by our clustering approach analysis. Our MD simulations at µs scales demonstrated the S-protein’s thermodynamics of the critical states at around 60°C, and the stable and denatured states for temperatures below and above this value, respectively.

Potentiometric Biosensors Based on Molecular-Imprinted Self-Assembled Monolayer Films for Rapid Detection of Influenza A Virus and SARS-CoV-2 Spike Protein.

Rapid, yet accurate and sensitive testing has been shown to be critical in the control of spreading pandemic diseases such as COVID-19. Current methods which are highly sensitive and can differentiate different strains are slow and cannot be conveniently applied at the point of care. Rapid tests, meanwhile, require a high titer and are not sufficiently sensitive to discriminate between strains. Here, we report a rapid and facile potentiometric detection method based on nanoscale, three-dimensional molecular imprints of analytes on a self-assembled monolayer (SAM), which can deliver analyte-specific detection of both whole virions and isolated proteins in microliter amounts of bodily fluids within minutes. The detection substrate with nanoscale inverse surface patterns of analytes formed by a SAM identifies a target analyte by recognizing its surface nano- and molecular structures, which can be monitored by temporal measurement of the change in substrate open-circuit potential. The sensor unambiguously detected and differentiated H1N1 and H3N2 influenza A virions as well as the spike proteins of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and Middle-East respiratory syndrome (MERS) coronavirus in human saliva with limits of detection reaching 200 PFU/mL and 100 pg/mL for the viral particles and spike proteins, respectively. The demonstrated speed and specificity of detection, combined with a low required sample volume, high sensitivity, ease of potentiometric measurement, and simple sample collection and preparation, suggest that the technique can be used as a highly effective point-of-care diagnostic platform for a fast, accurate, and specific detection of various viral pathogens and their variants.

Supervised machine learning approach to molecular dynamics forecast of SARS-CoV-2 spike glycoproteins at varying temperatures.

ABSTRACT: Molecular dynamics (MD) simulations are a widely used technique in modeling complex nanoscale interactions of atoms and molecules. These simulations can provide detailed insight into how molecules behave under certain environmental conditions. This work explores a machine learning (ML) solution to predicting long-term properties of SARS-CoV-2 spike glycoproteins (S-protein) through the analysis of its nanosecond backbone RMSD (root-mean-square deviation) MD simulation data at varying temperatures. The simulation data were denoised with fast Fourier transforms. The performance of the models was measured by evaluating their mean squared error (MSE) accuracy scores in recurrent forecasts for long-term predictions. The models evaluated include k-nearest neighbors (kNN) regression models, as well as GRU (gated recurrent unit) neural networks and LSTM (long short-term memory) autoencoder models. Results demonstrated that the kNN model achieved the greatest accuracy in forecasts with MSE scores over around 0.01 nm less than those of the GRU model and the LSTM autoencoder. Furthermore, it demonstrated that the kNN model accuracy increases with data size but can still forecast relatively well when trained on small amounts of data, having achieved MSE scores of around 0.02 nm when trained on 10,000 ns of simulation data. This study provides valuable information on the feasibility of accelerating the MD simulation process through training and predicting supervised ML models, which is particularly applicable in time-sensitive studies. GRAPHIC ABSTRACT: SARS-CoV-2 spike glycoprotein molecular dynamics simulation. Extraction and denoising of backbone RMSD data. Evaluation of k-nearest neighbors regression, GRU neural network, and LSTM autoencoder models in recurrent forecasting for long-term property predictions.